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INTRODUCTION 


Large amounts of sediment and other pollutants are carried annually 
In the rivers, lakes, estuaries, and coastal waters of the United States. 
These sediments and pollutants are major determinates of water quall':y. 
Many agencies are Investigating the potential of using remote sensing 
techniques to monitor various water quality parameters because of the 
ability of remote sensing to provide synoptic views over large areas. 

Specific data needs usually vary among different user o'gani- 
zations (Kuo, 1976b). Generally, the desired use of remote sensing 
data Is either identification or quantification of surface sediment 
and pollutants. This study Is concerned with data analysis procedures 
for quantification of water quality parameters that have already been 
Identified and are known to exist within the water body. Specifically, 
the study deals with the linear multiple-regression technique as a 
procedure for defining and calibrating data analysis algorithms for 
such instruments as spectrometers and multispectral scanners. The 
technique has been used by a number of authors (Johnson, 1975, 1977a, 
1977b, 1977c, 197?d, and Rogers, 1975, 1976) with apparent success. 
Unfortunately, results have not been completely satisfactory in that 
(1) analysis of subsets of data from the same experiment sometimes 
gives different correlations and algorithms (compare Johnson, 1975, 

1977a, 1977b); (2) repeated experiments over the same water body do 
not always allow quantification of the same water quality parameters 
(compare Rogers, 1975, 1976); (3) optimum results are not always 
obtained by a multiple regression equation with radiance as the 
independent variable (Ohlhorst, 1978 and Rogers, 1976); and (4) final 


mapped results do not always agree with results of other types of 
analysis (compare Bressette, 1978 and Johnson, 1977b). From these 
facts, It Is clear that a more complete understanding of the limita- 
tions, requirements, and precision of the linear multiple regression 
technique is required before It can be applied. by user agencies In an 
operational manner. 

In an effort to provide improved scientific understanding, an 
analytical and laboratory analysis of th« linear multiple regression 
technique has been conducted (Whitloc.^, 1977b) . That analysis deroon- 
b crated that the technique is fundamentally sounvi, and, in principle, 
should apply to many environmental situations in which both wecer and 
atmosphere contain linear and nonlinear optical effects. A number of 
limitations (optical, mathematical, and operational) were defined, 
however. 0n« problem with that analysis is that field results were 
not considered in that study. This investigation extends those results 
to include actual field data. An analytical investigation of the signal 
response equations is conducted, and results from "mixed-brew" laboratory 
experiments are presented. A study of measurement errors and ground- 
truth operations from several past field experiments is presented, 
and recommendations are made concerning future investigations from both 
aircraft and satellite instruments. 

SIGNAL- RESPONSE EQUATIONS 

Reflectance at sea level and upwelled radiance at altitude are 
related as follows (Miller, 1977): 

\ (X) - T^(X) [p^(X) (Lg^(X)) + L^^(X) + Lj.^(X)] + L^(X) 


( 1 ) 


where: 


L (X) > apparent upwelled radiance at altitude z at wavelength X. 
z 

T (X) > atmospheric transmission at wavelength X. 

P^(X) " inherent upwelled reflectance slightly above water surface 
at wavelength X. 

L (X) ■ upwelled radiance slightly above water surface from 100 
80 

percent diffuse reflector at wavelength X. 

L^j(X) ■ upwelled radiance from specular reflection of diffuse 
skylight at wavelength X. 

L (X) - upwelled radiance from specular reflection of sunlight 
t s 

(sun glitter) at wavelength X. 

L (X) ■ upwelled radiance from light scattered by the atmosphere 
8 

(path radiance) at wavelength X. 

For a remote sensing scene with a specific solar elevation angle, 
the inherent component of L^(X) being contributed by the water column 


is T^(X)[L^(X)] as follows: 


- P„(«a3„(x)) 

(2) 

T (X) IL (A)l . T (X)IP, (X)(L (X))] 

tt u & u so 

(3) 

where: 



L^(X) - inherent upwelled radiance slightly above water surface 
at wavelength X . 

A signal response model of the water may be assumed in which the 
remote sensing signal is expressed as the signal from the background 



water plus the signal change caused by some pollutant. The change 
caused by the pollutant may be expressed as a gradient constant times 
the pollutant concentration. The Inherent signal component for a 
simplified three-constituent water mixture may be assumed as follows: 

TglL^(X)] - A + BP^ ■»- EPg + (4-) 

where: 

A > ir erent upwelled radiance component from background water 
including loss due to atmospheric transmission. 

BP^ inherent upwelled radiance component caused by pollutant A. 

EPq > Inherent upwelled radiance component caused by pollutant B. 

SP^ ■> inherent upwelled radiance component caused by pollutant C. 

(assumed to vary nonlinearly as power Q) . 

B,E,S *• gradient constants includi.ng losses due to atmospheric 
transmission. 

P ,P ,P ■ concentrations of pollutants A, B, and C, respectively. 

A V w 

Upwelled radiance components from surface reflection and atmospheric 
path radiance effects car. be expressed as a value for the baseline 
atmosphere of a particular day plus a change caused by a variation 
in atmospheric pollutants over the scene as follows: 

where : 

1 ■ surface reflection and path radiance components for the 
baseline atmosphere over the scene. 



LX^ - change in surface reflection and path radiance conponents 

caused by atmospheric pollution (assumed to vary nonlinearly 
as power N). 

L ■ gradient constant. 

* concentration of atmospheric pollutant X^. 

In most remote sensing experiments, upwelled radiance is not measured 
at a specific wavelength, but instead an integrated average of radiance 
is measured over a range of wavelengths. To differentiate between 
wavelength specific and integrated-average values, the term Rad will 
he used to denote integrated-average values for apparent upwelled 
radiance at altitude Z over a range or band of wavelengths. Equation 
(1) may be rewritten as; 

Rad - A + BP, + EP_ + SpJ! + I + LX*J (6) 

ABC A 

If measurements are made with a remote sensing instrument with bands 
W,X,Y, and Z, then the equations for the measured values are: 

'‘“'Hi ■ *u * Va Vs ^ V? ^ 'w + Vl 

"“''x ■ *X * Va * Vs + v2 * 'x * vJ 

( 7 ) 

"“*'y * N Va Vc 'y Va 

■‘““z ■ *z Va + Vs * SjP? + i; * 


where subscripts W,X,Y and Z denote values over the same 
wavelength ranges as bands U,X,Y, and Z. 


Equation (7) is y. seirles ? four simultaneous equations which are 
linear in tour utkuowne ***** *■*** ®**^*'®®**^^<^*^ 

operations described In Api<ei.4lx C of reference by Whitlock (1977b) 
are folliwsd, equation (7) cau be solved to produce the following 
solutions for the values of concentration for pollutants A and B. 

- J + K,/Rad^) + + K^(Rad^) + K2<Rad2) (8) 

Pg » J’ + ^ 

The j* cons'.Ant*’ '.*ie a function of the A, B, E, S, I, and L 
constants .'f the '^^iciour hands, and the K, K' constants are a function 
of t'.f. B, L, S, and L constants. A key element in arriving at the 
above sr.lution is that the degree of nonlinearity in the P^ and 
contributions must be essentially constant over the wavelength range of 
bands W, X, Y, and Z or the changes must be small so that linear 
approximations are appropriate. 

Equations (8) and (9) represent an exact solution for two con- 
stituents which have linear radiance gradients in a remote sensing 
scene which (1) contains a water mixture with three constituents 
(one of which has a nonlinear radiance gradient), and (2) has nonlinear 
variations in surface reflection and path radiance contributions due to 
variation in some atmospheric constituent over the scene. While this 
model is somewhat simple, additional complications can be selectively 
added. Additional water pollutants (beyond pollutant C) and surface- 
atmospheric contribution terms can be added to equation (6) to compli- 
cate the model. As long as the degree of nonlinearity is constant 
over the wavelength range of Interest, a series of linear algebraic 


aquations will rasult and tha exact solution for and will 
be of tha same fom as aquations (8) and (9) except that ndditional 
bands of radiance values will be required. A key element of equations 
(8) and (9) is the values of P^ and P^ can be calculated without 
knowledge of each other or the values of P^ or X^. Knowledge of the 
A, B, E, S, I, and L constants of equation (7) is required, however, 
to compute values for the J, J', K, and K' constants before Rad 
values can be used to compute P^ and P^. 

In actual field experiments, values for the constants A, B, E, 

S, I, and L are seldom known. In fact, the number of water and 
atmospheric parameters which have a significant influence on the total 
upwelled radiance signal is usually unknown. As a result of this 
situation, empirical methods must be used to determine values for the J, 
J', K, and K' constants in eruatlons (8) and (9). Fortunately, 
equations (8) and (9) are of the same form as the statistician's 
traditional multiple-regression equation; 

p - 1 

P - J + I K. Rad + e (10) 

i - 1 ^ 

where : 

P dependent variable. 

Rad^ ■ independent variable. 

J,K^ * coefficients obtained by least squares fitting techniques. 

c ■ error. 

p - total number of J, K coefficients. ' 


In order to eveluate J, IL coefflcianta by empirical methods, 

certain euvironmental restrictions must apply, /dl data used in the 

fitting process must be from a situation with identical J, K 

coefficients. This means that values for the A, B, E, S, I, and L 

constants must be equal for each data point. This condition is often 

satisfied if one had a single remote sensing scene with multiple 

ground-truth points, and atmospheric transmission values do not vary 

by a large amount over the scene. (Surface reflections and path 

radiance are allowed to vary.) Values for the J, K coefficients 

generally will not be equal for each data point if one has the situation 

of multiple remote sensing scenes with a data point in each. In that 

case, each scene has a different solar elevation angle and atmospheric 

condition which causes different values for T (X), L (A), L_.(A), 

a so ro 

L (x), and L (A) in equation (1). This means that the A, B, E, S, 

f 8 A 

I, and L constants in equations (6) and (7) will be different for 
each data point. Since A, B, E, S, I, and L determine values for 
the J, K coefficients, the J, K coefficients will be different for 
each data point. In such a case, it is Impossible to obtain an 
accurate estimate of the J, K coefficients if radiance and ground- 
truth concentration values used in the least-squares multiple regression 
process are from different remote sensing scenes. This partly explains 
why linearized multiple regression with radiance as the independent 
variable failed to give a good data-reductlon algorithm in Ohlhorst 
(1978). In that experiment, the operational technique was multiple 
overpasses with a single ground-truth point in each scene. 


LEAST-SQUARES AND STATISTICAL METHODS 


EatiBmtlon of J, K CotfflclOTts - Th# ragrattlon cask Is to ascimata 
tha J, K coafflcianca in which Rad^ is assumad as Cha indapandanC 
variabla. In stany obsarvations, cha indapandanC variablas ara corralacad 
with aach ochar as wall as with cha dapandanc variabla which oakas 
results difficult CO incarpret (Snedacor» 1967). For cha ramota 
sensing situation, high corralacions betwaen the indapandanC variables 
(Rad^) should be expected if Cha pollutant of interest has a broad 
spectral signal over the wavelength range. 

As noted previously, least-squares procedures are used to estimate 
the J, K coefficients using a number of ground-truth points where 
radiance-constituent concentration data pairs are available. In 
performing the process on remote sensing data, three major assumptions 
are involved (Daniel, 1971): 

1. The correct form of the equation has been chosen (Rad^^ is 
linear with concentration for all bands involved). 

2. The data are representative of the whole range of 
environmental combinations in the remote sensing scene. 

3. The observations of the dependent variables (ground-truth 
concentration values) are uncorrelated and statistically 
independent. 

Three minor assumptions are: 

1. All observations of rhe dependent variable (concentration) 

2 

have the same (but unknown) variance, o . 


2. Tht distribution of uncontrollsd scror is norsksl. 

3. All indspsndent vsrisblss (Rsd^ vsluss) art known without 
arror. 

0ns problM is that naasurssMnts of ths Indspsndant variable (Rad^) do 

contain errors. Daniel (1971) indlcaras that errors in the Independent 

variable cause estiaates of the J, K coefficients to be biased. As 

a rule-of> thumb » it is recesnended that remote sensing experiments 

be designed such that the varicnce of radiance about mean values for 

2 

the ground truth locations ) l>e at least 10 times the variance 

2 ^ 

of data noise (o„ ). This rule-of>thumb may be referred to as Daniel's 
"i 

Criteria. In terms of standard deviation values: 

5 

°Rad, > 10.0 0 ^ (11) 

or 

Equation (12) states that ground-truth locations should be selected 

« 

within the remote sensing scene in such a manner that the standard 
deviation of the change in upwellcd radiance for the ground-truth points 
should be at least 3.16 times the standard deviation of data noise for 
the particular remote sensing instrument being used if least-squares 
procedures are to be used in the analysis of data. 

Measures of Precision - Before an experiment, it is not known how many 


bands (or in what combination^ will be required to separate the desired 


water-quality parameter from the total mix of factors which contribute 
to apparent upwelled radiance. The usual process is one of first 
calculating a regression equation for the best single band of radiance 
data and then succasslvely defining multiple-regression equations for 
the best two bands, three bands, etc. As additional bands are utilized, 
a number of statistical parameters may be used as indicators of the 
precision of each new multiple-regression equation. A number of factors 
must be considered when viewing these parameters to select an optimtuB 
multiple-regression equation. 

One popular statistical parameter used as a measure of precision 

is the correlation coefficient, r. The proportion of total variation 

2 

that is not explained by the regression equation is 1 - r . (An r 
value equal 0.9 means that 19 percent of the signal variation is not 
explained by the multiple-regression equation.) Draper (1966) stated 
that r is not a good measure of precision as the number of estimated 
coefficients approach the number of experimental observations. This 
implies that either the number of ground-truth observations should 
exceed the number of instrument bands by a wide margin or the number 
of bands included in the multiple-regression equation should be limited. 

The standard error, a, is a second measure of precision of the 
least-squares process for estimating J, K coefficients. The standard 
error is assumed to represent a value within which 68 percent of all 
errors are expected to fall if (1) there is an infinite number of 
observations, and (2) there is minimal error in the independent 
variables (Rad^) . Unfortunately, most remote sensing experiments 
have only a limited number of ground-truth observations and 
measurements do contain errors. 


The F-test Is a third method of evaluating the adequacy of the 

least-squares process. The calculated value of F must be greater 

than a critical value (F ) taken from F-dlstrlbutlon tables In order 

cr 

to be judged significant. If the multiple-regression equation Is to 

be used -for predictive purposes, calculated F should be at least A times 

the tabulated critical F value (Draper, 1966). One problem with this 

parameter Is that a level of confidence must be arbitrarily selected 

before F/F _ can be calculated, 
cr 

Daniel (1971) recommended the statistic, C^, as a measure of 

Che sum of the squared random errors. Given a multiple-regression 

equation with p estimated J, K coefficients, a low value of 

In combination with a C /p ratio < i.O is considered to indicate a 

P - 

good fit with negligible bias. The regression equation should then 

be useful for predictive purposes. The Cp/p ratio Is the only one of 

the above parameters which Is Indicative of bias in the fitted equation. 

For future remote sensing experiments. It is recommended that 

multiple- regression equations be computed for all combinations of bands 

for which upwelled radiance values are available. Values for r, o, 

F/F^^, Cp, and Cp/p should be computed for each equation. The prime 

basis for selecting an "optimum" regression equation should be 

C /p < 1.0 for minimum bias and F/F > 4.0 for predictive 
p — cr ~ 

utilization. Values for r should approach 1 and a should approach 
zero. Daniel's Criteria (eq. 12) should also be satisfied for all 
bands In the regression equation. Wlian all five conditions are met, 
the multiple regression equation with the minimum number of bands should 
be selected for calculation of pollutant concentration. 


LABORATORY VALIDATION EXPERIMENTS 


It is desirable to validate data analysis techniques with actual remote 
sensing data under controlled conditions. To achieve this result* tests 
were conducted with flltered-deionized tapwater in the Marine Upwelled 
Spectral Signature Laboratory at the NASA Langley Research Center. A 
sketch of the laboratory setup is shown in Figure 1. (More complete 
descriptions of the system and test procedures can be found in Whitlock 
(1977a and 1977b.) Only partial results from one series of tests will 
be presented in this paper for reasons of brevity. 

It was desirable to test the multiple regression technique with 
data from water mixtures which contained constituents with both linear 
and non-llnear radiance gradients. Single-constituent tests were con- 
ducted on a number of materials. From these data, it was concluded 
that both Ball Clay and Feldspar soils have near liner.r radiance 
gradients for concentrations between 4 ana 173 ppm. Rhodamine WT dye 
has a nonlinear gradient for concentrations between 17 and 1052 ppb. 

With this knowledge, a series of three-constituent tests were conducted 
with 25 different water mixtures. Table I shows the concentrations of 
Ball Clay, Feldspar, and Rhodamine WT that were present in the filtered- 
deionized tapwater for each test. Also shown are raiiance values for 
the following 5 wavelength bands: 


Band Number 

1 

2 

3 

4 


Wavelength Range 
(nm) 

240-500 

460-620 

540-700 

620-780 


Center Wavelength 
(nm) 

420 

540 

620 

700 


5 


700-860 


780 


The radiance values shown are in terns of relative units obtained by 

dividing power/bandwidth measurenents over the water by gray-card diffuse 

reflectance measurements of the input light source. 

A multiple regression analysis was performed for Ball Clay using 

12 of the 25 tests in Table I to simulate "ground- truth" values. The 

12 "ground- truth" values were tests 1,‘ 3, 5, 6, 8, 10, 13, 15, 18, 20, 

21, and 23. The standard deviation of the change in upwelled radiance 
• ' 
was next calculated for each band for this "ground-truth" data set. 

Values are compared with data noise from the laboratory measurements 

as follows: 


Band Number 

°Rad 

(Wliitlock, 1977b) 


1 

0.1372 

0.0343 

4.00 

2 

0.1125 

0.0343 

3.23 

3 

0.1414 

0.0343 

4.12 

4 

0.1431 

0.0343 

4.17 

5 

0.0992 

0.0343 

2.89 


These figures indicate that Daniel's Criteria (eq. (12)) is satisfied 

for all bands except number 5. Since the signal-to-nolse standard 

deviation ratio for band 5 is only slightly below the value of 3.16, 

it was .lecided not to exclude that band from the regression analysis 

in this particular case. Regression equations for all combinations 

of bands were next calculated for Ball Clay concentration. Estimated 

values of the J, K coefficients and various statistical measures of 

precision are shown in Table II. Review of Table II Indicates that 

the lowest value of total squared error (C ) is obtained for band 

P 

combinations 2, 3, and 4. The regression equation for this combination 


is: 


(13) 


- -A.l + 234.4 (Rad 2 ) - 613.7 (Rad^) + 918.0 (Rad^) 

statistical estimates of precision for this equation are: 

r - 0.99 
0 ■ 6.8 ppm 

C - 3.0 
P 

C /p ■ 0.8 
P 

Since values for r« o, and C^/p are acceptable, it is assumed 

that equation (13) has good predictive capability. To test this 
assumption, the radiance values for bands 2, 3, and 4 from Table I 
were applied to equation (13). Values for Ball Clay concentration 
(P^) were calculated and are compared with actual values in Figure 2. 
Shaded symbols denote the 12 "ground-truth" points used in the least- 
squares fit, and the open symbols represent calculations for the 
remaining 13 independent test points. (Some open symbols are hidden 
under the shaded points.) Since all points fall within + 3.9o of the 
true value, it was concluded that equation (13) has good predictive 
capability and that linear multiple-regression analysis procedures 
have the potential for quantification of constituents with linear 
radiance gradients in water mixtures which also contain nonlinear 
constituents. A satisfactory linear regression equation was also found 
for the concentration of Feldspar soil in these same mixtures (Whitlock, 
1977b). 


REVIEW OF FIELD DATA 


Analytical analysis and laboratory test casas have been used to 
perform a limited validation of linearised multiple-regression analysis 
for quantification of marine constituents. No matter how many controlled 
\ tests are conducted, final validation of the technique must come through 

use of field experiments. As noted previously, results from past 
experiments have not been completely satisfactory. This section examines 
the data and operational conditions of past tests in an attempt to 
define problem areas which require increased attention for future 
experiments. 

Instrument Noise - Unfortunately, remote sensing data always contain 
error and random noise. Most operational instruments contain onboard 
calibration lamps and black bodies to minimize radiance error. Noise 
has been somewhat more difficult to eliminate because it is caused by 
many components in the instrument system. After an experiment, instru- 
ment noise can usually be evaluated by calculating the standard deviation 
of calibration lamp and/or black body count values and applying the 
appropriate calibration constants to convert to radiance. If onboard 
calibration sources are not available, flat or near-constant radiance 
portions of the remote sensing scene may be examined. Both of these 
procedures have been used to examine noise in the data for the field 
experiments described in Johnson (1977b, 1977c, 1977d) and Ohlhorst 
(1978). Figure 3 presents estimates of noise standard deviation values 
for each band for which data were obtained. These data indicate that a 
wide range of noise values have been observed in previous experiments. 


The problem of data noise has been widely recognized, and pixel 
averaging is often employed as a method of reducing the error. To 
illustrate the effects of this process, calibration lamp values from 
the experiment (Johnson, 1977b) were averaged in a manner to simulate 
pixel averaging. New values for noise standard deviation were calculated 
and are compared with single-pixel values in Figure 4. Large reductions 
in instrument noise may be achieved by averaging even small size (4 by 4 
or less) pixel arrays. Use of array sizes larger than 7 by 7 produced 
only small reductions in noise for the data set (Johnson, 1977b). 

Pixel averaging f'.ffectively increases the size of the ground 
resolution element that is observed by the remote sensing Instrument. 

If the array size becomes too large, average radiance values of the 
enlarged pixels may not be representative of ground-truth values because 
of spatial variability in pollutant concentration. As pixels are 
averaged, noise error will be reduced but new error is introduced because 
of scene dynamics or hydraulic changes. Carried to extreme, hydraulic 
features that depict smaller-scale pollutant transport may be erased 
by the averaging process, "'f averaging is necessary, then the following 
criteria should be used to establish credibility of the enlarged pixel 
radiance values. 


^°Rad ^nxn ^°N^nxn 
Pt 


(14) 


where : 

^**Rad ^^nxn 
pt 


standard deviation of enlarged pixel about average 
radiance value (calculated from individual pixel 
radiances within n by n array surrounding each ground 
truth point) . 




When [Oq.j ]_ is much larger than [a„] , ona haa an indication 

REa_. nxn N nxn 

Pt 

that scene dynamics have introduced large errors into the process that 
offset the noise teductlon benefits gained by pixel averaging. 

Ground Truth Placement - Previous discussion has indicated that ground- 

truth points should be located within the remote sensing scene in such 

/ 

a manner that Daniel's Criteria (eq. 12) is satisfied for all bands in 
the regression equation. To review this aspect of previous experiments, 

f 

I 

values of have been calculated using radiance values over each 

set of ground-truth points for the field experiments of Johnson (1977b, 
1977c, 1977d) and Ohlhorst (1978). Pixel averaging had been employed 
to eliminate signal changes caused by the ground-truth boat, time 
differences, and location uncertainty as well as to reduce noise in each 
data set. New estimates of o„ have been made from calibration lamp 
or background water data to simulate the pixel averaging process for 
each experiment. Values of have been computet end are com- 

pared with Daniel's Criteria in figure 5. It appears that ground-truth 
placement was adequate such that Daniel's Criteria was satisfied after 
pixel averaging for most experiments. Approximate calculations indicate 
that not a single one of these experiments would satisfy Daniel's 
Criteria on a 1 by 1 pixel array basis, however. Ground-truth locations 
for future remote sensing experiments should be selected in such a 
manner that the largest possible value to minimize the 

requirement for pixel averaging. Use of past remote sensing Images 
of the area to estimate relative radiance differences between proposed 
locations is probably the best method of selecting a distribution with 


maximum standard deviation. 


Physical Consistency - An unwritten assuaption In the miltlple regression 
process Is that all data used In the correlation are "good" without 
physical deficiencies. When remotely monitoring water-quality parameters, 
this generally means that a number of criteria should be satisfied 
concerning both environmental conditions and physical operations of 
the experiment. Only limited discussion justifying physical con- 
sistency of the data is available in Johnson (1977b, 1977c, 1977d), Rogers 
(1975, 1976) and Ohlhorst (1978). For chat reason, the remainder of this 
paper will discuss some of the practical problems experienced by the 
authors 

Relative to environmental criteria, two desirable conditions 

are: 

1. Water depth greater than the Secchi depth. 

2. Constant vertical concentration gradient within the remote 
sensing penetration depth (the depth above which 90 percent of 
the upwelled radiance originates.) (Gordon, 1975). 

It is sometimes assumed that the maximum remote sensing penetration 
depth is the same as the Secchi depth. Calculations in McCluncy, (1974) 
as well as unpublished data from the James River in Virginia Indicate 
that maximum remote sensing penetration depth may be on the order of 
20 to 50 percent of Secchi depth, depending on absorption and scattering 
characteristics of the mixture. Knowledge of the maximum remote sensing 
penetration depth is required so that water samples are not obtained 
below the zone th"t is causing the remotely sensed signal. One method 
to estimate the maximum remote sensing penetration depth is that of 
lowering a flat black plate with less than jL percent diffuse reflectance 
into the water and noting its depth of disappearance (McCluney, 1974). 


In terns of opsrstlonsl critsris, s major problem is the tlM 
lapse between overpass of the remote sensor and collection of individual 
water samples. Time lapse can cause the ground- truth data to be 
hydraulically inappropriate if there are significant wind or tidal 
Influences on the water body being observed. It may be possible to 
correct ground- truth data to account for small time lapses (Kuo, 1976a), 
but such procedures have not yet been widely demonstrated in field 
experiments. In principle, all data should be synchronous with the 
remote sensor overpass, but that has not been achieved in past experi- 
ments, as shown below: 


Reference 

Geographical 

Area 

Number of 
Ground Truth 
Points in 
F egression 
Analysis 

Maximum Time 
Lapse Between 
Water Sampling 
and Remote 
Sensor Overpass 

Johnson, 1977b 

James River in Virginia 

21 

2.0 hrs 

Johnson, 1977c 

New Jersey Coast 

22 

3.0 hrs 

Johnson, 1977d 

New York Bight 

10 

0.3 hr 

Rogers, 1975 

Saginaw Bay in Mlchingan 

27 

8.0 hrs 

Rogers, 1976 

Saginaw Bay in Michigan 

16 

8.0 hrs 

Ohlhorst, 1978 

Delaware Shelf 

7 

0.5 hr 


Some experiments have a small time lapse, but in other cases it has been 
assumed that constituents In a water pixel remain constant for 2 to 8 
hours. Considering that algae tend to migrate depthwise with changing 
light intensity (solar elevation angle) , large time lapses should be 
justified with quantitative water sample data even when wind and tidal 
effects are such that flow conditions nre stable. 


In the operational area, another major problem la often 
encountered. Large-scale experiments with multiple vessels often require 
either some water sample analysis onboard ship or laboratory analysis 
of the same parameter by different organizations. Consistency of 
laboratory results between different laboratories is a longstanding 
problem, and variations of certain parameters may bo larger chan that 
of the remote sensing scene. It is recommended that future experi- 
ments maintain consistent handling and laboratory analysis procedures. 

PRACTICAL APPLICATIOHS 

Axi Investigation of the linear multiple-regression technique with 
remotely sensed radiance as the Independent variable has been conducted. 
Slgn.il response equations have been analyzed, and results from "mixed 
brew" laboratory tests are presented. Results from these studies 
Indicate that the technique is fundamentally sound and should apply 
in many environmental situations in which both the water and atmosphere 
contain linear and nonlinear optical effects. Conditions which limit 
application of the tecinique have also been discussed. A review of 
previous field experiments has served to emphasize additional limita- 
tions which must be considered in future experiments. A summary of 
recommended conditions for use of the technique is given in Figure 6. 

From this listing, It is clear that the linearized multiple-regression 
analysis should never be applied blindly to a set of data without back- 
ground knowledge concerning atmospheric variability over the scene, 
the constituent of interest, hydraulics of the water body, typical 
Secchi and maximum remote sensing penetration depth values, and iMasure- 
ment uncertainties from various sources. The technique has strong 


theoretical foundation and careful application should jrlsld useful 
restilts. It la particularly appropriate for uae In eaell regions 
for single- tlaM experlaenca to validate either ■atheaatlcal ur hydraulic 
aodels of pollutant transport and diffusion. Present cost for 
obtal**lng mltlple ground-truth points within a single rwote sensing 
scene Halts Its usefulness for aany routine aonitorlng alssiona, 
towever. 

CONCLUSIOHS 

The study described herein Is part of a continuing effort to 
define data-reduction techniques at d their appropriate application so 
that increased benefits can be derived from both aircraft and satellite 
renote sensing data. The goal of the present study was a nore complete 
understanding of limitations, requirements, and precision of the linear 
multiple-regression technique with radiance as the independent variable. 
Environmental a id optical ptiysics conditions have been defined for 
which an exact solution to the signal response equations is of the sas» 
mathematical form as the statistician's traditional multiple-regression 
equation. In such a case, use of linearized multiple regression is 
merely an empirical correlation to obtain coefficients for the exact 
solution to the signal-response equations. Additional analytical 
Investigations are desirable to more completely define atMspheric 
limitations and to consider the problem of bottom reflection in 
optically shallow waters. 

One problem with the use of the regression technique is that the 
independent variables (upwelled radiances) contain errors and are often 


corrs.'atsd with Mch oth«r. Thlt raquirta eonaidaratlon of • nui^or of 
•catittical paraaotora whan parfondng tha ragraaalon analyaia. 

Ravi aw of paat flald axparlaanta Indlcacaa that data noiaa vaa of 
auch aagnltuda that data aaoothini waa ntquirad bafora Danial'a Critarla 
could ba aatlafiad for tha laaat-aqiiaraa aultlpla-ratraaaion proeaaa. 
laprovtd ailactlon of groimd-truth locationa to laaxlaizo varianca la 
racomandad to alnlalzo data aaoethlng raquiraoanta and tha phyalcal 
arrora aaaociatad with that proeaaa. 

Tiaa lapaa batwaan raaota aanaor ovarpaaa and watar aa^la 
collaetion appaara to hava baan a problaa in paat axparlaanta. Econoaie 
conaidaration will alwaya raault in only a liadtad nuabar of watar aaaplaa 
obcalnad aynchronoualy with tha rtaota aanaor ovarpaaa. Additional 
atudiaa to davalop and daaonacraca ttchnlquaa for corractlng non~ 
synchronous data for raaota sensing use are daslrabla. 
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Table I. Concentration and radiance data for three constituent tests 
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Table II. Estimates of J, K coefficients and measures of precision 
for Ball Clay sediment in three-constituent mixtures. 
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Figure 4. Effects of pixel avei'aglng for data In Johnson (1977b) 
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Figure 5. Signal to noise standard deviation ratio for several field experlaencs. 
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Figure 6. Rec«-ended condltlone for use of linear ■ultlple-regreeeloo technique 
with reliance aa the Independent variable. 


